Stacking-fault energies for Ag, Cu, and Ni from empirical tight-binding potentials 



O 

o 

(N 



(N 



R. MeyeiQ and L. J. LewisQ 

Departement de physique et Groupe de recherche en physique et technologie des couches minces (GCM), 
Universite de Montreal, CP. 6128 Succursale Centre-ville, Montreal (Quebec) H3C 3J7, Canada 

(Dated: July 12, 2002) 

The intrinsic stacking-fault energies and free energies for Ag, Cu, and Ni are derived from 
molecular-dynamics simulations using the empirical tight-binding potentials of Cleri and Rosato 
[Phys. Rev. B 48, 22 (1993)]. While the results show significant deviations from experimental 
data, the general trend between the elements remains correct. This allows to use the potentials for 
qualitative comparisons between metals with high and low stacking-fault energies. Moreover, the 
effect of stacking faults on the local vibrational properties near the fault is examined. It turns out 
that the stacking fault has the strongest effect on modes in the center of the transverse peak and 
its effect is localized in a region of approximately eight monolayers around the defect. 

PACS numbers: 61.72.Nn,63.20.Dj,02.70.Ns 
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I. INTRODUCTION 

The stacking-fault energy is an important property of 
a material since it determines _to a high degree its de- 
formation and failure behaviorE Therefore, for the com- 
puter modeling of processes involving plastic deformation 
and/or failure, a correct representation of the stacking- 
fault energy by the employed model is desirable. How- 
ever, many models have been constructed without consid- 
eration of stacking-fault energies and the predicted va 
deviate often considerably from experimental results! 
A precise knowledge of the stacking-fault energies pre- 
dicted by empirical models is therefore of high interest 
in order to correctly interpret results of computer simu- 
lations. 

The .-empirical tight-binding potentials of Cleri and 
RosatcCI are extensively used in atomistic simulations 
since they give a reasonable description of many struc- 
tural and thermal properties of fee metals. In this work 
we study the intrinsic stacking-fault energies and free 
energies predicted by these potentials for Ag, Cu, and 
Ni and compare them with the results of similar models 
for Cu and Ni. The stacking-fault energies were calcu- 
lated at T = 300 K from molecular-dynamics simulations 
and at T = K with the help of relaxation calculations 
(molecular-statics). The derivation of the vibrational 
density of states (VDOS) from the simulations made it 
also possible for us to calculate the stacking-fault free 
energies within the quasiharmonic approximation and to 
study the local effect of stacking-faults in fee metals. 

The aim of this paper is to provide the stacking-fault 
energies for Ag, Cu, and Ni modeled by the potentials 
of Cleri and Rosato. We expect that these values will 
contribute to proper discussions of results based on these 
potentials. Although in the case of Ag and Cu the val- 
ues derived from the simulations are significantly lower 
than the experimental results, the trend between all three 
elements is well reproduced. It has to be noted that 
the occurrence of low stacking-fault energies is a general 
problem of central-force potentials. One reason for this 
is that without the consideration of angular forces sig- 



nificant contributions to the stacking-fault enerey arise 
from the third and higher neighbor shells, only.EI A re- 
lated explanation is given by the nearest neighbor d-band 
tight-binding model where the energy difference between 
the fee and hep structure is dominated.-by the fifth and 
sixth moment of the density of states p whereas many 
empirical potentials, including those of Cleri and Rosato, 
are based on the second moment approximation. Better 
agreements of computationally simple empirical poten- 
tials with experimental values for the stacking-fault en- 
ergy require therefore the explicit consideration of the 
stacking-fault energy during the construction of the po- 
tentials. 



II. COMPUTATIONAL METHODS 
A. General 

As stated in the introduction, the stacking-fault ener- 
gies and free energies reported in this article have been 
derived from molecular-dynamics simulations employing 
the potentials of Cleri and Rosato.El In these simulations 
the equations of motions have been integrated with, the 
help of the velocity form of the Verlet-algorithmil us- 
ing a time-step of 2 fs. Isothermic-isochoric (NVT en- 
semble) and isobaric-isoenthalpic (NpH ensemble) condi- 
tions were achieved where necessary by the application of 
the Nose-Hooverpthermostatcl method and the Parrincllo- 
Rahman schemeE3 (restricted to orthogonal simulation 
boxes), respectively. Periodic boundary conditions were 
applied in all cases. 

Calculations of the relaxed structures at T = K were 
done in two steps. First, the systems were cooled down 
to a temperature of T — 0.1K. Afterwards, the final 
configurations were obtained with the help of a steepest- 
decent method. 

For each of the three elements a stacking-fault configu- 
ration and a reference configuration were prepared. The 
stacking-fault configurations were composed as a stack of 
110 triangular lattice layers building a nearly cubic sys- 
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tern of N st = 1 018 160 atoms with one stacking-fault; the 
reference configurations on the other hand are regular fee 
systems A^f cc = 1 000 188 atoms (63 3 cubic cells). 

At first glance, configurations with more than a million 
atoms may seem to be exaggerated for the simple calcu- 
lation of stacking-fault energies. However, these system 
sizes guarantee the presence of a sufficiently large num- 
ber of vibrational modes for the calculations of free en- 
ergies in the quasiharmonic approximation. In addition, 
the large system sizes lead to large distances between the 
periodic copies of the stacking-fault, thereby minimizing 
possible self-interaction effects. 

Before the calculation of the stacking-fault energies 
could begin, the equilibrium dimensions of all configu- 
rations had to be determined at T — and 300 K. In 
the case of the reference configurations these were de- 
rived from the T = K lattice constants given in Ref. || 
and the average system volumes obtained from isobaric- 
isoenthalpic simulations at T = 300 K and zero pressure. 

For the determination of the equilibrium sizes of the 
stacking-fault configurations, it had to be taken into ac- 
count that only volume relaxations perpendicular to the 
plane of the stacking fault are physically relevant. Relax- 
ations in the directions parallel to the fault-plane are arti- 
facts of the finite system size and diminish if the number 
of atomic layers parallel to the fault is increased. For this 
reason the in-plane dimensions of the simulation boxes of 
the stacking-fault configurations were derived from the 
lattice constants obtained for the perfect systems and 
only the out-of-plane dimensions had to be explicitly de- 
termined. At T = 300 K these missing dimensions were 
again obtained with the help of zero-pressure isobaric- 
isoenthalpic simulations in which only the out-of plane 
direction of the simulation box was allowed to fluctuate. 
Similarly the volume relaxations at T = K were derived 
from relaxations of the stacking-fault configurations with 
fixed in-plane dimensions of the simulation box. 



energy or enthalpy, this procedure may lead to a small 
offset in the temperature which in turn affects the value 
of the potential energy. While these discrepancies usually 
are very small and can be neglected in most cases, they 
can have an dramatic impact on the values of 7 since 
the potential energy difference on the right hand side 
of Eq. (Q) is much smaller than the potential energies 
themselves and can be comparable to — or even smaller 
than — the error in the potential energies induced by 
temperature deviations. In order to avoid this problem 
we derived the potential energies from simulations in the 
canonical (NVT) ensemble. 



C. Calculation of stacking-fault free energies 

The exact calculation of free energies from molecular- 
dynamics simulations is a difficult task. In—this work 
we employ the quasiharmonic approximation^ to derive 
the free energies F of the stacking-fault and reference 
configurations. Within this approach the free energy is 
given by 



F = E (V)+W / g(w) 



hui 1 , 
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where j3 is the inverse temperature (k^T) , g(u>) the 
normalized VDOS and Eq(V) denotes the ground-state 
energy of the system for the equilibrium volume at tem- 
perature T. In contrast to most molecular-dynamics sim- 
ulations, Eq. (||) takes quantum effects fully into account. 
However, since we are mainly interested in the values 
of the stacking-fault energies in a classical computer- 
simulation context, it is more appropriate to calculate 
the free energy in the classical limit. Taking this limit 
and neglecting the purely quantum zero-point fluctua- 
tions, Eq. (0) becomes 



B. Calculation of stacking-fault energies 



F cl = E (V) + 3Np- 



g(ui) \og((3huj)diLj. (3) 



After the dimensions of the configurations had been de- 
termined, the stacking-fault energies 7 of the three metals 
could be obtained. For this the simple formula 
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was used, where A is the area of the stacking-fault plane, 
N s f the number of atoms in the stacking-fault configura- 
tion and e s f, ef cc the potential energies per atom in the 
respective configuration. The kinetic energy does not 
contribute to the stacking-fault energy since its per par- 
ticle value does not depend on the configuration. 

Particular attention has to be paid to the system tem- 
perature. The usual method to attain the desired tem- 
perature in molecular-dynamics simulations is to rescale 
the particle velocities repeatedly until the system is equi- 
librated. In simulations run under conditions of constant 



Given Eq. ^) the task of calculating the free energy 
of the configurations by molecular-dynamics simulations 
reduces to the determination of the vibrational density 
of states g(u>) and the ground-state energy Eq(V). The 
latter quantity was obtained by a relaxation of the sys- 
tems to T — K using the simulation box dimensions of 
T = 300 K. In order to derive the vibrational density of 
states g(u)), simulations in the NVE ensemble were per- 
formed at T = 300 K. From the trajectories of the parti- 
cles during these simulations the velocity-autocorrelation 
function (v(t)v(O)) was derived which is related to the 
normalized VDOS byB 



(v(t)v(O)) 
_ x (v(O)v(O)} 1 



t dt. 



(4) 



In order to reduce statistical errors, the density of states 
g{uj) was averaged over ten or more simulation runs of 
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5 000 steps, each. From the density of states obtained 
this way the free energies of all six configurations could 
be calculated with the help of Eq. (||) . These values were 
then used to calculate the intrinsic stacking-fault free en- 
ergies 7 of the three elements by replacing the potential 
energies per particle in Eq. (0) by the corresponding free 
energies per particle. Since at T = OK energy and free 
energy are identical these calculations had to be done 
only at T = 300 K. 



III. RESULTS 
A. Stacking-fault Energies 

In Table | we present the stacking-fault energies and 
free energies obtained from molecular-dynamics simula- 
tions as described in the preceding section. From this 
table it can be seen that the general agreement with the 
experimental results is only qualitative, similar to compa- 
rable models. The calculated stacking-fault energies (and 
free energies) for Cu and Ag are much smaller than their 
experimental counterparts; the calculated stacking-fault 
free energy of Ag at room temperature is even negative. 
For Ni the situation is not clear due to the large variation 
of the experimental data. However, the trend among the 
three elements is well reproduced by the model. This 
means that while one has to be careful with the interpre- 
tation of quantitative results, the potentials of Cleri and 
Rosato are well suited for comparisons between materials 
with high, low and very low stacking-fault energies. 

In accordance with the experimental data,E3 the results 
in Table [j| show a decrease of the stacking-fault energies 
and free energies with increasing temperature. In the 
case of the stacking-fault energies, the decrease is caused 
by anharmonic effects only, since the specific heat of a 
purely harmonic, classical system does not depend on 
structural details. In contrast to this, the decrease of 
the stacking-fault free energies in Table | contains har- 



TABLE I: Experimental (e) and theoretical (t) values of 
stacking-fault energies for Ag, Cu, and Ni (in mj/m 2 ). 7 and 
7 are the stacking-fault energies and free energies obtained in 
this work; OJ and VC denote tk values given in Ref. M-fe*i 
the models of Oh and Johnsonta and Voter and ChcnljO 
The experimental resuiH* are taken from the compilations of 
GallagheiO and Murr.O 
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FIG. 1: a) Calculated normalized VDOS g{ cc (ui) for crystalline 
fee Cu. b) Difference Ag(cu) — g B f(t-u) — gi cc (u) between the 
VDOS g S f(u) of the system with stacking-fault and gf cc (uj). 



monic and anharmonic contributions. However, due to 
the nature of the quasiharmonic approximation, the an- 
harmonic contributions are treated only in an approxi- 
mate fashion via the volume dependence of Eq and the 
temperature dependence of the VDOS. This approxima- 
tion gets in the way if one wants to compare the stacking- 
fault energies and free energies in order to determine the 
role of entropic contributions to the stacking-fault en- 
ergy. Particularly striking is the large discrepancy be- 
tween the stacking-fault energy and free energy of Ni at 
room temperature. Since our results indicate that the 
entropy difference favors the stacking-fault configuration, 
the stacking-fault free energy should be the lower quan- 
tity, as it is the case for Cu and Ag. This shows that, at 
least in the case of Ni, there are important anharmonic 
contributions to the stacking-fault free energy which are 
not accounted for by the quasiharmonic approximation. 



B. Vibrational Density of States 

Figure [l] shows the VDOS calculated for the fee struc- 
ture as well as the difference of the VDOS between the 
stacking-fault and the perfect fee configuration. It should 
be noted that the absolute value of this difference is phys- 
ically irrelevant since it depends on the system size and 
goes to zero in the limit of an infinitely thick system. 
However, the shape of this function is system size inde- 
pendent and reveals the energy regions where the VDOS 
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FIG. 2: Normalized local VDOS of the atoms in the first, 
second, and fourth layer away from a stacking fault in Cu 
compared to the VDOS in the perfect crystal. 

is affected by the stacking fault. From the lower panel of 
Fig. [I] it can be seen that the strongest effect of the stack- 
ing fault appears in the transverse peak where modes are 
transferred from the center of the peak to (mostly) lower 
energies. A similar but smaller transfer takes place in 
the longitudinal peak. It is these transfers which give 
rise to the entropy difference between the configurations, 
which in turn accounts for the temperature dependent 
change of the stacking-fault free energy within the quasi- 
harmonic approximation. A similar behavior is found for 
Ag and Ni with a slightly more pronounced transfer of 
longitudinal modes in the case of Ni. 

In order to see the spatial extent of the influence of 
stacking faults on the vibrational properties of the crys- 
tal lattice we have calculated the local VDOS for the 
first atomic layers near a stacking fault in Cu. The re- 
sults of this calculation, which are presented in Fig. ||, 
show that atoms in the first and second layer near the 
fault account for most of the change of the total VDOS 
in Fig. [l], whereas the VDOS in the fourth layer differs 
only slightly from the bulk VDOS. Atoms in the third 
layer (whose VDOS is not given in order to avoid over- 
crowding of the figure) have an intermediate character. 
While their VDOS is still clearly different from the VDOS 
of the crystalline bulk, it is already much closer to this 
than to the VDOS of the first two layers. 

IV. SUMMARY AND CONCLUSIONS 

The stacking-fault energies predicted by the empiri- 
cal tight-binding potentials of Cleri and Rosato reflect 



the experimental values of Ag, Cu, and Ni qualitatively. 
Similar to other empirical models they underestimate the 
values of Ag and Cu. In contrast to this and in rough 
agreement with experimental data, rather high stacking- 
fault energies are found for Ni. Since the experimental 
data for Ni show strong variations it is impossible to say 
whether the potential for this metal overestimates the 
stacking-fault energy or not. In any case, the potentials 
can be used for qualitative comparisons of materials with 
high, low, and very low stacking-fault energies. 

For all three metals we find a decrease of the stacking- 
fault energy with increasing temperature of the same or- 
der of magnitude as it is experimentally observed. til This 
decrease has to be considered in the construction or eval- 
uation of models which explicitly take the stacking-fault 
energies into account. 

The discrepancy between the calculated room tem- 
perature stacking-fault energy and free energy for Ni 
shows that important anharmonic effects are missed by 
the quasiharmonic approximation. For this reason one 
should not put too much emphasis on the negative value 
found for the room temperature stacking-fault free en- 
ergy of Ag. What remains is the question why the an- 
harmonic effects are so much stronger for Ni. A possible 
explanation might be given by the relative importance of 
the many-body term vs. the pair-potential term in the 
model. An inspection of the parameters of the poten- 
tials in Ref. |^ shows indeed that the pair-potential term 
has a lower weight and a shorter range for Ni than for 
Cu or Ag. Thus, the increased weight of the many-body 
term in the case of Ni might account for the increased 
occurrence of anharmonic effects in this metal. 

Finally, the analysis of the local effect of stacking- 
faults on the vibrational properties in fee metals reveals 
a strong transfer of modes to lower energies in the trans- 
verse peak and a similar but less pronounced transfer 
in the longitudinal peak of the VDOS. Inspection of the 
layer-resolved local VDOS shows that these effects are 
localized in a region of six to eight monolayers around 
the defect. 
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